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RESEARCH MEMORANDUM 


A METHOD FOR EVALUATING THE LOADS AND CONTROLLABILITY 

ASPECTS OF THE PITCH-UP PROBLEM 

By Melvin Sadoff, Frederick H. Matteson, 
and C. Dewey Havill 


SUMMARY 


A procedure is described for estimating the range of peak airplane 
load factors and maneuvering tail loads likely to be experienced in 
pitch-up maneuvers. The method assumes a realistic evaluation maneuver 
which partially integrates airplane and pilot response. 

Results of computations, in which it is assimied that this evaluation 
maneuver is used on an example swept-wing airplane at 35^000 and 15^200 
feet, indicated that though the load factors and maneuvering tail loads 
were not critical in pitch-up maneuvers at 35^000 feet, they were likely 
to exceed design levels at 15,200 feet. It was shown, however, that for 
corrective-control rates of ^5^ per second or higher, the airplane load 
factors would not exceed the design value by more than 10 percent. It 
was indicated that it would be desirable to restrict the maximum avail- 
able corrective-control rate to some optimum value which compromises the 
high rates required to minimize overshoot load factor with the low rates 
desirable for low maneuvering tail loads. 

A tentative criterion, based on the ratio of the destabilizing moment 
at the time of corrective-control application to the corrective-control 
moment per unit stick deflection available to the pilot, appears promising 
for predicting controllability of pitch-up. Preliminary information on 
two swept-wing fighter airplanes indicated that ratios of 1 to 20 and 1 
to 100 were associated with an uncontrollable pitch-up and a relatively 
controllable pitch-up, respectively. 


INTRODUCTION 


One of the important problems experienced with swept-wing airplanes 
is the undesirable pitch-up tendency associated with nonlinear pitching- 
moment curves. The pitch-up is considered undesirable in two main 
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respects. Under certain conditions^ the design wing and horizontal-tail 
loads may be exceeded inadvertently in pitch-up maneuvers. Furthermore, 
the occurrence of pitch-up limits controlled maneuvering, in many cases, 
to load factors below the pitch-up boundary; in others , it results in a 
significant reduction in controllability. 

Several flight investigations (refs. 1 to k) and an analytical study 
(ref. 5) provide some information of general interest in connection with 
the loads and controllability aspects of the pitch-up problem. Refer- 
ence 1 presents some experimental evidence of the possibility of encounter 
ing large wing and tail loads in pitch-up maneuvers. In references 2 to 
4, pilot opinion of the pitch-up characteristics and controllability of 
two 35^ swept-wing airplanes with various wing and tail modifications is 
presented. The analytical investigation of reference 5 assesses various 
pitching-moment irregularities in terms of the abruptness of the airplane 
response during pitch-up to a more or less arbitrary control input; how- 
ever, no specific attempt was made to consider quantitatively the load 
factors and tail loads or the relative controllability that may be 
expected in specified pitch-up maneuvers. 

The purpose of the present paper is to outline a procedure, based on 
an assumed realistic control input (or evaluation maneuver) , which may be 
used to assess primarily the loads aspects of the pitch-up problem. Also, 
the possibility of predicting the degree of controllability of a speci- 
fied pitch-up tendency is briefly discussed. Computations are made for 
an example swept-wing fighter airplane to illustrate the use of the method 


NOTATION 


— + M A , 

airplane damping coefficient, — — - — 1/sec 

mV 1-w- 


Ms 




Cq control-deflection coefficient, ^ l/sec‘ 


Z5 

control-rate coefficient, — l/sec 

mV 


Cm 


airplane pitching-moment coefficient about airplane center of 
airplane pitching moment 
qSc 


gravity. 


Cl 


airplane lift coefficient. 


qS 
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airplane normal-force coefficient, 


wing mean aerodynamic chords ft 


g acceleration of gravity^ 1 g = 32.2 ft/sec^ 


pressure altitude, ft 


airplane pitching moment of inertia, slug-ft^ 


K parameter denoting damping ratio of airplane to that of 

horizontal tail 


-Mi-v / '■ 

airplane spring constant, l/sec^ 


Z-b distance from airplane center of gravity to aerodynamic center 

of horizontal tail, ft 


airplane lift , lb 


I>|. horizontal- tail lift, lb 


AMq intercept of a particular linear segment of the pitching-moment 
curve on the ordinate axis, Aa = 0 


m 


W 

airplane mass, — , slugs 


N 


airplane normal force, lb 


n airplane normal load factor , ^ 


dynamic pressure, Ibs/sq ft 
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wing area, sq ft 


S-j; horizontal-tail area, sq ft 


variable introduced in Laplace transform 


time, sec 


V airplane velocity, ft/sec 


W 


airplane weight, lb 


a airplane angle of attack, deg or radians 


flight-path angle, radians 


elevator angle, deg or radians 


Sstick control-stick angle (fore and aft), radians from neutral 


when preceding symbol denotes increment from steady-state 
condition 


downwash angle, deg or radians 


% 

horizontal-tail efficiency factor, 


angle of pitch, radians 


p mass density of air, slugs/cu ft 

dCr , 

Cl airplane lift-curve slope, 1/ radian 
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(Cl«) 


Wt 


horizontal-tail lift-curve slope, , l/radlan 






Cmc 


Ma 


M, 


a 


Mj= 


M. 


Mi 


M« 


-•a 


^5 


l/ radian 

d6e 


l/radian 

da 


l/radian 

dSe 


(CniQ_) q.SC, ft-lb/radian 


%(M)’ 


ft-lb 


radian/sec 


qSc, ft-lb/radian 

Sstick ft-lb/radlan 




pVStIt' 

t ft-lb 


2 > radian/sec 

ft-lb 


radian/sec 

iVradian 


qsj, Ib/radian 

<x, 7 } ®e>\ 

5, n i 

a, 7 , 0 equivalent notation for and 

^ ^ ^ dt^^ dt2" dt^ 


. ^ ^ „ da ^7 d0 , dn 

equivalent notation for — ^ and ^ 


d^a d^7 


d^0 
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Subscripts 


t 


horizontal tail 


e 


corresponding to a specified value of pitching acceleration 


i 


initial conditions in ith time interval 


max 


maximum value 


th 


threshold 


D 


duration of pitching-acceleration stimulus 


DESCRIPTION OF EXAMPLE AIRPLANE 


The example airplane used in the present computations is a jet-powered 
swept-wing fighter type. A photograph of the airplane is presented in 
figure 1, and the physical characteristics and a two-view drawing of the 
airplane are given in table I and figure 2, respectively. 


A realistic evaluation of the loads that might occur during a pitch-up 
requires the establishment of a rational pitch-up evaluation maneuver. 

In general, a pitch-up maneuver may be expected to consist of three dis- 
tinct parts: (l) application of elevator control at a constant rate until 

the pilot detects pitch-up; (2) continued application of control at the 
initial rate for an interval depending on the pilot's reaction time; and 
( 3 ) application of corrective control to arrest the pitch-up at a rate 
depending on several factors including pilot experience with pitch-up, 
intensity of pitch-up, and proximity to the design load factor. 

It is necessary in order to arrive at a rational evaluation maneuver 
to identify and determine the level of a response quantity which the pilot 
associates with the onset of pitch-up. It is also necessary to determine 
a reasonable average pilot reaction time. Some limited results in refer- 
ences 2 and 3 indicated that an appropriate response quantity might be 
pitching acceleration, since the pilots appeared to associate different 
levels of pitch-up intensity with the magnitude of pitching acceleration 
developed during pitch-up. In view of this, tests were conducted on a 
modified Link trainer, to ascertain whether there existed a minimum value 


DESCRIPTION OF METHOD 


Evaluation Maneuver 
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of pitching acceleration which correlated consistently with the pilots* 
initial appreciation of pitch-up. The results of these tests ^ which are 
described in more detail in Appendix A^ indicated a fairly consistent 
correlation of a pitching acceleration threshold A0 of about O.I 5 with 
the pilots* initial perception of pitch-up. The tests also indicated a 
mean reaction time of about O .3 second. For the present analysis, an 
additional 0.1 second (or a total of 0.^ second) was applied to account 
for the time required to accelerate the control surface from rest to a 
constant rate. 

With the concept of a pitching-acceleration threshold established , 
and a reasonable reaction time, and with the pertinent aerodynamic and 
geometric data known, it was then possible to define parts (l) and (2) of 
the evaluation maneuver. Since it was not possible to predict the exact 
control rate used by a pilot to arrest a specified pitch-up tendency, a 
range of corrective-control rates v/as selected for part ( 3 ) of the evalua- 
tion maneuver to illustrate the effect of this variable on the load fac- 
tors and tail loads that may be experienced in pitch-up maneuvers . 

The control inputs, established by the above procedure for the example 
airplane, are presented in figure 3 for a Mach number of O.9O (the speed 
at which the pitch-up is most severe for the example airplane) and for two 
altitudes - 35,000 feet and about 15,000 feet. The upper altitude was 
chosen to correspond to the altitude at which most of the flight tests were 
performed on the example airplane and the lower altitude was selected to 
illustrate the loads aspects of the pitch-up problem at low altitude . In 
the present example, this altitude corresponds to the altitude at which 
the pitch-up flight region (lower boundary defined herein as the angle of 
attack for neutral stick-fixed stability) was just penetrated in a 6 g 
maneuver. The several initial control rates (fig. 3)^ which correspond to 
initial values of A of 0.2, 0.5, and 1.0 g per second, were selected to 
cover a reasonable range of entry rates from the relatively gradual maneu- 
vers used by Ames pilots in research tests to the more abrupt maneuvers 
that are likely to be used in training or combat. 


Computational Procedure 


The equations of motion and the pertinent aerodynamic and geometric 
data used in the computations are presented in Appendix B. Airplane 
responses associated with the elevator motions shown in figure 3 were 
obtained with a Reeves Electronic Analog Computer. 

In the event that an analog computer is not readily available, or 
where there are only a few configurations to check, a computational pro- 
cedure using the Laplace transform method is also described in Appendix B. 
A sample set of computations is presented and the results are compared 
with the solutions obtained from the REAC . 
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APPLICATION OF METHOD TO EXAMPLE 
SWEPT-WING FIGHTER AIRPLANE 

Loads 


Computed response quantities for the example airplane at 35^000 feet 
are presented in figure 4. Incremental angle-of-attack and pitching- 
acceleration variations are shown in figure 4(a), while computed time 
histories of airplane load factor and incremental maneuvering horizontal- 
tail load are given in figure 4(b). The dashed lines in figures 4(a) and 
4(b) marked ^'lirnit of ACjj^(a) and ACj^(a) curves*’ refer to the maximum 
values of Aa and n for which experimental pitching-moment and lift data 
used in the computations were available. The two vertical ticks on the 
Aa and n responses indicate, respectively, the time at which the pilot 
would first perceive pitch-up (A0 = O.I 5 ) and the time at which corrective 
control was initiated. The results in figure 4 indicate that for 
corrective-control rates of 20^ per second or less, recovery does not 
generally occur within the angle-of-attack range for which the computa- 
tions are available, and the overshoot^ values of Aa and n are relatively 
large, exceeding 9^ ^*5^ respectively. For the higher corrective- 

control rates, recovery occurs at incremental angles of attack of less 
than 8^ and at load factors less than 4, and the overshoot in angle of 
attack and load factor is about 3-5^ ^-nd 1.5^ respectively. There does 
not appear to be any consistent effect of entry rate n on the peak air- 
plane responses during pitch-up (fig. 4). However, it may be pointed out 
that for a constant corrective-control rate, the overshoot values of angle 
of attack and load factor increase appreciably with an increase in entry 
rate. It should also be noted that the values of Aa and n at the 
pitching-acceleration threshold decrease appreciably with an increase in 
entry rate. The results in figure 4 also show little variation of the 
peak negative pitching acceleration or maximum positive incremental maneu- 
vering tail^load with entry rate n. In the present example, the peak 

values of G and AL^.. for corrective-control rates above 20^ per second 

G 

are limited by the maximum down-elevator deflection available. 

Figure 5 presents computed response quantities for the example air- 
plane in pitch-up maneuvers at 15^200 feet. Incremental angle-of-attack 
and pitching-acceleration time histories are presented in figure 5(3*)^ Q'^d 
normal load factor and incremental maneuvering horizontal-tail load varia- 
tions are shown in figure 5 (l). Comparison of these results with those of 
figure 4 indicates that for the same values of corrective-control rate and 
entry rate, the overshoot in angle of attack at the lower altitude is only 
about 60 percent of that at 35^000 feet. In the present case, this results 

^In the present analysis, overshoot is defined as the difference 
between the peak values of Aa and n and the values existing at the time 
the threshold value of pitching acceleration was attained. 
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in overshoot load factors at 15,200 feet about 50 percent greater than 
those at 35 >000 feet. It may be noted that the peak load factors reached 
in the assumed pitch-up maneuvers at 15,200 feet exceed the design value 
appreciably for the lower corrective-control rates. However, for 
corrective-control application at a rate of ^5^ per second or higher, the 
value of n^g^ would not exceed the design level by more than 10 percent. 
For the same corrective-control rates, the peak negative pitching accelera- 
tions and the corresponding maximimi incremental tail loads attained in 
pitch-up maneuvers at 15,200 feet (fig. 5) approximately 50 to 100 per- 
cent greater than the peak values reached at 35 >000 feet (fig. 4). The 
results in figure 5(^) show, for a constant corrective-control rate, rela- 
tively small difference in due to changes in entry rate n, 

®max 

while a substantial increase in tail load is associated with an increase 
in corrective-control rate. These and the foregoing results suggest that 
it may be desirable to limit the maximum available control rate to some 
optimum value which compromises a reduction in overshoot load factor with 
a decrease in the maneuvering tail load. 

Summary plots , which may be useful in selecting an optimum control 
rate, are presented in figures 6 and 7 for altitudes of 35,000 and 15,200 
feet, respectively. Values of overshoot load factor and incremental maneu- 
vering horizontal- tail load are plotted as a function of corrective-control 
rate for various values of n. The results at 35,000 feet (fig. 6) indi- 
cate that for corrective-control rates above ^5^ per second, only a small 
further reduction in overshoot load factor occurs, while practically no 
change occurs in the values of ALt.. . At 15,200 feet (fig. 7), for 

^max 

corrective-control rates greater than about ^5^ P^^ second and up to the 
maximum rate considered of 75^ P^^ second, the overshoot load factor is 
further reduced by only approximately 0.2 (3 percent of the design value), 
while the values of increase from roiighly 90 percent to 120 

^max 

percent of the design tail load for the example airplane. It appears, on 
the basis of these results, that it might be desirable to limit the maxi- 
miorn available control rate to about 45^ per second, since a further 
increase in rate results in an appreciable increase in maneuvering tail 
load without materially reducing the overshoot load factor. 


Controllability 


In a previous flight study of the pitch-up problem on the example 
airplane (ref. 2), it was indicated that pilot opinion of the pitch-up 
appeared related to the level of peak positive pitching acceleration 
experienced during the pitch-up maneuvers. Althoiigh the magnitude of the 
peak pitching-acceleration response may describe the relative controlla- 
bility of a pitch-up tendency on a given airplane, where the effectiveness 
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of several modifications on the pitch-up is being evaluated (as in ref. 2 , 
for example) , this response quantity alone is not sufficient to indicate 
the relative controllability of pitch-up on several different airplanes 
with different pitching moments of inertia and longitudinal control effec- 
tiveness. In a preliminary assessment of the problem, it was decided that 
a more general controllability criterion might comprise all three of these 
quantities, that is, pitching acceleration, moment of inertia, and control 
effectiveness , and would relate the unstable pitch— up moment ^y^max 
the time of corrective-control application to the corrective-control 
moment per unit stick deflection available to the pilot ^Sstick* 


A plot presenting the variation of the controllability parameter 

. with entry rate h is shown in fig\ire 8 for the example 
s"bick 

airplane at 35,000 and 15,200 feet. It may be noted that an increase in 
value of this petrameter implies a decrease in controllability of the asso- 
ciated pitch-up tendency. The results in figure 8 indicate a significant 
reduction in controllability with an increase in entry rate n for both 
35,000 and 15,200 feet. It also appears from the results in figure 8 that 
the pitch-up for the example airplane is more easily controlled at 15,200 
feet than at 35,000 feet, since the corrective -control moment per \mit 
stick deflection increases more rapidly than the destabilizing 

moment lySmax with a decrease in altitude.^ 


Flight tests at 0.90 Mach number at 35,000 feet, where entry rates 
n of 0.2 g per second to 0.5 g per second were used to enter the pitch-up 
region, have indicated that the pitch-up tendency on the \Himodified exam- 
ple airplane was relatively uncontrollable. From the upper curve in 
figure 8, it may be noted that the computed values of controllability 
factor corresponding to this uncontrollable pitch-up vary between about 
0.05 and 0.065. In order to provide some information on the magnitude of 
the parameter ly^max/^bg-j^^Qj^ corresponding to an airplane which has a 

relatively mild pitch-up and which is considered fairly controllable, com- 

putations were also made for the 3$^ swept-wing airplane described in 

2It sho^Ild be recognized, however, that since the controllability 

parameter I,,^' „„/Mc ^ , is roughly an inverse measure of the ability 
^ J max/ Ostick ^ /-u n i 

of the pilot to reduce a given destabilizirg moment (.hence angular accel- 
eration) to zero , the improved controllability in the present case refers 
only to the ability of the pilot to control airplane attitude, that is, 
angle of attack, angle of pitch, etc. This was touched upon previously 
in the discussion where it was noted the overshoot in angle of attack at 
15,200 feet was only about 60 percent of that at 35^000 feet, (if it is 
desired to define controllability as the ability of the pilot to control 


load factor, then the proposed controllability parameter should be multi- 
plied by appropriate values of dynamic pressure and lift-curve slope. 

The controllability parameter would then be greater at the lower altitude, 
i.e., a reduction in controllability. This is in agreement with the 
larger overshoot in load factor at the lower altitude previously noted.) 
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reference 4. This airplane differs from the example airplane mainly in 
that its pitching moment of inertia and longitudinal control effectiveness 
are greater than those of the example airplane by factors of 1-1/2 and 4, 
respectively. The computed results for this airplane are shown in fig- 
\ire 9 where they are compared with the data for the example airplane at 
35,000 feet. It may he seen that the values of the controllability param- 
eter for the reference airplane are only about 20 percent of those for 
the example airplane over the rarge of entry rates n covered in the com- 
putations. 

It is recognized that the proposed controllability criterion should 
be checked for a number of airplane configurations for which pilots' 
opinions of the pitch-up characteristics are available. A study of this 
kind, directed toward prediction of pitch-up characteristics, is currently 
under way. 


CONCLUSIONS 

A method is described for predicting the rajige of airplane load fac- 
tors and horizontal-tail loads that may be experienced in pitch-up maneu- 
vers. The method is based on a realistic evaluation maneuver wherein it 
is assumed that the pilot applies nose-up longitudinal control at a con- 
stant rate until 0.^ second after the pitching acceleration has increased 

0.15 radian per second per second above the steady-state value, at which 
time nose-down corrective control is applied at various constant rates. 
Application of the procedure to an example swept -wing fighter airplane has 
led to the following conclusions: 

1. At 35^000 feet where the pitch-up region was entered at load 
factors well under the design value, the results of the analysis indi- 
cated that the airplane load factors and the incremental maneuvering 
horizontal-tail loads likely to be experienced in pitch-up maneuvers were 
not critical. The load factors would generally be restricted to values 
below the design level, either by the application of adequate corrective 
control or by the stall, while the tail loads were limited, in general, 
by the maximum down-elevator control available. 

2. At 15,200 feet where the pitch-up region was entered in a maneu- 
ver between 5 S 6 g, the analysis indicated both the load factors and 
tail loads likely to be encountered in pitch-up maneuvers would exceed 
design values. However, for corrective -control application at the rate 
of 45 ^ per second or higher, the peak load factors would not exceed the 
design level by more than 10 percent. 

3 . Based on the results of the present analysis, it appears desir- 
able to restrict the maxim-um control rate on airplanes which experience 
pitch-up to a value based on a compromise between desirable reductions 
in overshoot load factors and maneuvering tail loads. For the example 
airplane considered in the present analysis, an optimum rate of 45^ per 
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second is indicated, since greater rates would result in greater tail 
loads without materially effecting a further reduction in overshoot load 
factor . 


k. A computed parameter, based on the ratio of the unstable pitch-up 
moment at the time of corrective-control application to the corrective- 
control moment per unit stick deflection available to the pilot, appears 
promising for predicting controllability of pitch-up • Preliminary informa- 
tion available on two swept-wing fighter airplanes indicated that for a 
pitch-up considered uncontrollable by the pilots, the ratio ^y^max/^Sg^^^^j^ 

is of the order of 1 to 20, while for a relatively controllable pitch-up, 
the ratio is about 1 to 100. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Apr. 6, 1955- 
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APPENDIX A 
LINK TRAINER TESTS 


These tests were designed to determine empirically the values of 
pilots* reaction time and pitching-acceleration threshold needed to con- 
struct the pitch-up evaluation maneuver used in this report. The problem 
is analogous to that described in several research projects directed 
toward defining the transfer functions of a human operator. The available 
literature and results are summarized in references 6 to 9* Although 
these references contained some results on reaction time to a pitching- 
acceleration stimulus and on threshold values of pitching acceleration^ 
these data were not, in general, considered applicable to the present 
analysis since they were not obtained in an appropriate environment. 

In order to keep the environment as close to flight as practicable 
the tests were made in a Link trainer modified so that an outside operator 
could introduce arbitrary pitching motions which were independent of those 
associated with operation of the control stick. The trainer, figure 10, 
was equipped with instruments to record stick position, pitching velocity 
and acceleration, and time. The procedure was to measure the reaction 
time separately and then to deduce the pitching-acceleration threshold 
from the type of run shown in figure 11. As it was not possible to measure 
the reaction time to a vestibular stimulus using the modified Link trainer, 
the reaction time to a visual stimulus was used because, in the present 
case, the proper environment was considered more important than the type 
of stimulus used. 

The technique for measuring reaction time was to have the pilot make 
a gradual pull-out. At various random positions of the Link the operator 
would turn on an indicator light « The pilot was instructed to push the 
stick forward immediately upon seeing the light. The time interval between 
the light going on and the start of control application was considered the 
reaction timco This procedure was also used with the pilot distracted by 
having to hold a constant heading, with and without rough air. The test 
results from 20 to 30 runs each on three pilots and three engineers are 
summarized in table II. The mean reaction time is about 0.3 second which 
is in good agreement with the various data cited in reference 6. 

As indicated on figure 11 the procedure in the runs to determine 
pitching-acceleration threshold was to have the pilot make a gradual pull- 
out; then at a random time and using various angular velocity and accelera- 
tion rates, the outside operator would induce a pitch-up. The pilot was 
instructed to recover by pushing forward on the stick as soon as he per- 
ceived a pitch-up. As previously shown the threshold pitching accelera- 
tion was determined as the value existing a fraction of a second, equal 
to the subject *s mean visual reaction time, before the application of 
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corrective control. (See fig. 11 «) In order to establish whether visual 
perception or angular position or velocity was a factor^ runs were made 
with the pilot looking at the instrument panel with the gyro horizon 
covered and uncovered^ with the pilot looking out through a window in the 
canopy^ and with the pilot *s eyes closed. Tests were again run with the 
pilot being distracted by the necessity for maintaining a constant heading 
with and without simulated rough air. No appreciable effect was noted due 
to any of these variables. 

The results from 20 to 30 runs on each of three pilots and one engi- 
neer indicated a mean threshold of 0.12 radian per second as shown in 
table III. Since both the reaction time and threshold are statistical 
quantities there was a fairly large scatter indicated by the standard 
deviations. As shown in figure 12 the value of 0.12 radian per second per 
second is in general agreement with results from reference 9 which were 
presented as a function of the time duration of the stimulus . It may be 
noted that the pitching-acceleration threshold tends to decrease with an 
increase in the duration of the stimulus up to about 8 seconds . A limited 
series of flight tests (4 runs) on an F-84-F airplane, which were made to 
check the Link data, indicated a mean value of 0.l8 radian per second per 
second. The value used to construct the model pitch-up evaluation maneu- 
ver was selected as 0.15 radian per second per second. 
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APPENDIX B 

METHODS OF COMPUTATION 
REAC Method 


If a constant speed maneuver is assumed and if some of the higher 
order derivatives are neglected, the longitudinal equations of motion used 
may be written as 


-mV/ = Z(a) + Zg^ASe 

(Bl) 

lyO = M(a) + M^d + + Mg^ASg 

(B2) 


where the nonlinear functions Z(a) and M(a) were obtained from the upper 
two solid curves in figure 13 ^ and the other pertinent information was 
obtained from tables I and IV. It may be noted that either for the sake 
of simplicity or because insufficient information was available to define 
the variations, the values of M^, M^, and M^^ were assumed constant 

over the angle-of -attack range. Solutions were then obtained of Aa(t) , 
0(t), and /(t) for the longitudinal control inputs shown in figure 3- 


Laplace Transform Method 


In this method the nonlinear pitching-moment and lift curves are 
divided into several linear segments approximating the original curves, 
as illustrated in figure 13* Computations are then made for several time 
intervals where a new interval is dictated either by a change in slope of 
the linearized pitching-moment curve (fig. 13 ) or by the application of 
corrective control. The longitudinal equations of motion used in this 
case are 


-mV7 = ZoAa + (B3) 

Iy§‘ = MoAo. + + M^e + Ms^ASg (b4) 

where and are appropriate to the particular linear segment 

of pitching-moment curve under consideration. (See fig. 13*) The values 
of Zg^, Mgg, and were assumed constant over the entire range of angle 
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of attack. Since it is assumed that 7 = 0 - a, equations (B3) and (bU) 
may be reduced to the equivalent second-order equation 


a + ha + kAx = COA60 + + {/Mq/Ij) ^ 


(B5) 


where ASg = beC't - ti) + A6ej^ and (AMQ/ly)^ is proportional to the inter- 
cept of the linear pitching-moment segment under consideration in the ith 
time interval on the ordinate axis (Aa = O) . ^The general Laplace trans- 
formation of equation (B5) neglecting the Cj&e term^ which is generally 
small, and for a step (of magnitude plus a ramp elevator motion may 

be expressed as 


s^Zkt,(s) + sbAa(s) + kAa(s) - (s + b)Aa^ - 


. CQbe 

= “^i" 


CoASe^ + (AMo/ly)i 

s 


from which 


(s + b)Aaj_ + dj_ ^ 

^ s2 4. bs + k s(s2 + bs 4- k) s2(s2 -f bs 4- k) 

(b6) 

and 

% 

(s 4- b)Aous 4- d^s 

~ s2 4. bs + k s2 + bs + k s(s2 + bs + k) 

(B7) 

The b, k, Aa^, and d-i values are appropriate to the particular linear 
segment of pitching-moment curve being considered. The term CQAbe^ 
is proportional to the change in pitching moment due to elevator deflection 
from the start of the maneuver (Aa = 0; t = O) to the beginning of the I'C 
time interval. The airplane response in the time plane (inverse tr^s- 
formation of eqs. ^6) and (B?)) may be readily evaluated by Heavisides 
partial fractions expansion. (See ref. 10.) 

To illustrate the use of the Laplace transform method for computing 
airplane response in pitch-up maneuvers, a sample set of computations is 
presented for the example airplane at 35^000 feet using ^the longitudinal 
control input shown in figure ll^- for an initial value n of 0.2 g per 
second. This control input is slightly different from the corresponding 
one used in the REAC analysis because of a slightly different initial 
pitching-moment slope. (See fig. 13.) The pertinent basic data used are 
given in table IV and in figure 13* 
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For region I, which corresponds to the linear pitching-moment segment 
between values of Aa of 0 and O.OLL5 (fig. 13)^ 


Aa( s ) = 


Co^e 


0.195 


32(32 + bs + k) s^(s^ + 2.2s + 28.6) 


(B8) 


and 


a(s) 

The inverse transformations 
Aa(t) = -0.00052 

d( t ) = 0 . 00683 


= 0-195 

s(s^ + 2.2s + 28.6) 

of equations (B8) and (B9) are 
+ 0.00683t + an oscillatory tern 
+ an oscillatory term 


(B9) 


(BIO) 

(Bll) 


The oscillatory contributions in equations (BIO) and (Bll) are generally 
small and may usually be omitted . 

For region II (O.O628 > Aa > 0.0445) ^ 


_ 0.0445sQ + 0.1048sg + 0 . 490 s + O-I 95 
s2(s^ + 2.2s + 10.7) 


(B12) 


or 


Aa(t) = 0.0422 + 0.0l82t + (O.OO7I cos 3.085t - O.OO89 sin 3-085t) 

3.005 

(BI3) 

and 

d(t) = 0.0182 + (-0.0353 cos 3-o85t - 0.0121 sin 3-085t) 

(B14) 

For region III (O.O768 > Aa > O.O628), 


Aa( s ) 


0.0629 s 3 + 0.i293s2 - 0.452s + 0.195 
s2(s2 + 1.7s - 8.0) 


(B15) 
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or 


Aa(t) = 0.0514 - 0.0244t + 0.0154e^’^^^ - 0.0038e“3 (B16) 


and 

a(t) = -0.0244 + 0 . 03256 ^ + O.Ol45e"®*®^‘'^ (BI 7 ) 

Similarly for region IV ^ which corresponds to the linear se^ent of 
pitching-moment curve between values of Aa of O.O 768 and 0.21 (limit 
of ACm(a) and AC]^(a) curves)^ 


, , 0.076733 + o.i849s^ - 0.264s + 0.195 

Aa( s ) 

s 2 (s^ + 1 . 7 s - 6.6) 

or 

Aa(t) = 0.0323 - 0.02961 + 0.0447e^'®®'^ + 0.0002e"® 

and 

d(t) = - 0.0296 + 0 . 0832 e^-®®^ - 0 . 0007 e“ 3 * 56 t 


(BI 8 ) 


(BI9) 


(B20) 


At the point in time where corrective control is applied (fig. l4) 
the terms in the numerator of equation (B18) will change (since the 
initial conditions and the elevator ramp rate have changed) , but the 
denominator will remain the same (corrective control applied in region IV 
in present example) so that for a corrective-control rate of 10° per 
second. 


or 


Aa( s) 


0.1l48s^ + 0.34 o6s^ - 0.l86s - 1.62 

s^(s^ + 1.7s - 6.6) 


(B21) 


Aa(t) = 0.091 + 0.2455t - 0.00294e^*®®^ + 0 .0267e"° (B22) 


and 


d(t) = 0.2455 - 0 . 005476^*®®* - 0.095e“®*®®^ (B23) 

For corrective-control rates of 20°, 45 ^ and 75, per second, it is 
only necessary to substitute respective values of Cq 50 of -3.24, -7*3^ 
and - 12.17 for -1.62 in the numerator of equation (B21) . The results 
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of "ths compu.'ts.'bions for this case are presented in table V and in fig- 
ure iL, where they are compared with the corresponding solutions from the 
REAC. The agreement shown is fairly good, although it is indicated that 
linearizing the pitching-moment curve results in a slightly lower peak 
positive pitching acceleration (and therefore lower overshoot in angle of 
attack) than the values obtained from the REAC solutions. 
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TABLE I.- DESCRIPTION OF EXAMPLE AIRPLANE 


Wing 

Total wing area, sq ft 287.90 

Span, ft 37.12 

Aspect ratio 1-.79 

Taper ratio O.51 

Mean aerodynamic chord, ft 8.O8 

Dihedral angle, deg 3.O 

Sweepback of 0.25-chord line, deg 35-23 

Geometric twist, deg 2.0 

Root airfoil section (normal to 0.25-chord line) . . . NACA 0012-64 

(modified) 

Tip airfoil section (normal to 0.25-chord line). . . . NACA 0011-64 

(modified) 

Horizontal tail 

Total area, sq ft 3^»99 

Span, ft 12.75 

Aspect ratio 4.65 

Taper ratio 0.45 

Dihedral angle, deg 10.0 

Mean aerodynamic chord, ft 2.89 

Sweepback of 0.25-chord line, deg ..... 3^*59 

Airfoil section (parallel to center line) NACA 0010-64 

Maximum stabilizer deflection, deg 1 up, 10 down 

Elevator 

Area, sq ft 10. 13 

Span, each, ft 5-77 

Maximum deflection, deg 35 up, 17 *5 down 

Boost hydraulic 

Horizontal-tail length, ft 18.25 

Airplane weight, lb 12,400 

Airplane mass, slugs 385 

Airplane pitching moment of inertia, slug-ft^ 17^480 

Center-of -gravity location, percent c 22.5 
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TABLE II.- SUMMARY OF VISUAL REACTION TIME TESTS 


1 

1 

Pilot 

Engineer 


Item 

A 

B 

c 

A 

B 

c 

Number runs 

33 

19 

22 

53 

17 

19 

Mean reaction 
time, sec 

0 , 3 i ^3 

0.247 

0.246 

0,325 

0.283 

0.31 

Standard devi- 
ation, sec 

0.085 

0.045 

0.048 

0.125 

0.036 

0.117 


table III.- SUMMARY OF PITCHING-ACCELERATION-THRESHO^ TESTS 



Pilot 

Engineer 

Item 

A 

B 

C 

B 

Number runs 

20 

30 

19 

22 

Mean pitching- 
acceleration threshold , 
radians/ sec^ 

0.093 

0,143 

0.120 

0.124 

Standard deviation 
radians/ sec^ 

0.046 

0.090 

0.072 

0.055 


TABLE IV.- PERTINENT AERODYNAMIC DATA USED IN ANALYSIS 

Airplane lift curve, ACuCa) See fig. 13 

Airplane moment curve, ACm(a) 

Horizontal-tail lift-curve slope, per radian. . 3*0 

Elevator moment effectiveness, Cmg , per radian .... -0.24 

Downwash factor, Se/Sa fiS* 13 

Ratio of horizontal tail to wing dynamic pressure, 

/ Assumed 1.0 

Ht/‘l ! 

Damping ratio of airplane to that of horizontal tail, ^ i or 

^ .... 35,000 

Pressure altitude, hp, ft 87s 

Airplane velocity, V, ft/sec 

Mach number . O.9O 

Mass density of air, p, slugs/cu ft 0.00073® 
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TABLE V.- SAMPLE COMPUTATIONS 


® 

© 

© 

© 

© 

© 

© 

© 

© 

t 

Aa 

• 

a 

a 

a 

. y 

“^d 

mV 

•• 

e 

© +© 

n 

ALtg 

^t 

Control 

rate 

0 

0 

0 

0 

0 

0 

0 

0 

Entry 
- 1 . 21 ^/sec 

6.6 

.0445 

.00683 

0 

.007 

.007 

2.29 

-7 

7.75 

.0629 

.0223 

-.008 

.023 

.015 

2.68 

-l 4 

8.5 

.0767 

.0545 

.146 

.024 

^.170 

2,94 

-163 

8.55 

.1148 

.1454 

CVJ 

00 

.063 

.387 

3.46 

-370 

8.75 

.153 

.191 

1 — 1 

• 

.083 

.258 

3.93 

-250 

Recovery 
10^/ sec 

9.15 

*^.232 

.218 

0 

.095 

.095 

4.90 

-90 

8.75 

.148 

.161 

-.131 

.070 

-.061 

3.86 

60 

Recovery 

9.15 

.181 

-.040 

-.87 

-.02 

-.890 

4.27 

850 

9.35 

.152 

-.261 

-1.33 

-.11 

- 1.440 

3.92 

1380 

9 . 5 ^+ 

.0785 

-.564 

-1.92 

-.24" 

-2.160 

2.51 

2070 

20^/ sec 

9.565 

.0631 

-.613 

- 2.04 

-.27 

®-2.3l 

2.47 

2220 

9.585 

.0534 

-.645 

-1.63 

-.67 

-2.30 

2.38 

2200 

8.75 

.143 

.087 

-.850 

.038 

-.812 

3.80 

777 

Recovery 

45 °/sec 

8.95 

. 13 i^ 

-.192 

-1.900 

-.080 

-1.980 

3.70 

1890 

9.114 

.074 

-.576 

- 2.840 

-.250 

f-3.090 

2.88 

2960 

8.75 

.138 

-.003 

-1.711 

0 

-1.711 

3.74 

l64o 

Recovery 

75 °/sec 

8.85 

.126 

-.221 

- 2.640 

- .096 

" -2.736 

3.59 

2620 

8.889 

.114 

-.333 

-2.997 

-.144 

^- 3.141 

3.46 

3015 


^Obtained by solving equation (B 5 ). 


b 


-■a 


appropriate to linear pitching-moment segment under consideration. 


^Threshold ^ 
^Aa > limit of 
^Peak negative 
^Peak negative 


of 0.15 reached. 

ACj^(a) and ACm(a) curves. 

^ attained . 

^ reached at maximum available down-elevator deflection. 
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Figure 1.- The example swept-wing airplane. 
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(a) 35^000 feet. 

Figure 3 .- Longitudinal-control inputs used in analysis. 
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0 t for n = 0.2 ^20 22 24 

t, sec 

(b) 15,200 feet. 

Figure 3»- Concluded. 
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A a , deg 
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f ACm(a) and ACmCo) curves— 7 
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(a) Aa and 0. 

Figure 4.- Computed response quantities at 35^000 feet. 
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(b) n and AL-t-^ . 
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Figure 4,- Concluded. 
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( a ) Aa and 6 . 

Figure 5*- Computed response quantities at 15^200 feet. 
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(b) n and AL-j^* 
Figure 5 »~ Concluded. 
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^ *^overshoot 



AL.v 


lb 



5p , deg /sec 

^recovery ’ 


Figure 6.- Variation of overshoot load factor and peak incremental 

horizontal- tail load with elevator-control rate used in recovery at 
35,000 feet. 
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Figure 7«>“ Variation of overshoot load factor and peak incremental 

horizontal-tail load with elevator-control rate used during recovery 
at 15^200 feet. 
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n 

Figure 8.- Variation of controllability paramater Iy^max/%gtick 
entry rate n at 35 >000 feet and 15 >200 feet. 
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®max 

“stick 

radians 



Figure 9*~ Comparison between the controllability factors for the 
elevator- controlled example airplane with those for a similar 
stabilizer-controlled airplane at 35^000 feet. 
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Figure 11.- Procedure used to determine pitching-acceleration threshold. 
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Figure 12.- Pitching-acceleration thresholds from present tests and from 

tests of reference 9* 
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Figure 13.- Pertinent basic data used in analysis; pressure altitude, 

35,000 feet. 
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Figure l4.- Comparison of the angle-of-attack and pitching-acceleration 
responses obtained from the REAC and by the Laplace transform method . 
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